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Abstract. Grimmett's random-orientation percolation is formulated 
as follows. The square lattice is used to generate an oriented graph such 
that each edge is oriented rightwards (resp. upwards) with probability 
p and leftwards (resp. downwards) otherwise. We consider a variation 
of Grimmett's model proposed by Hegarty, in which edges are oriented 
away from the origin with probability p, and towards it with probability 
1 — p, which implies rotational instead of translational symmetry. We 
show that both models could be considered as special cases of random- 
oriented percolation in the NE-quadrant, provided that the critical value 
for the latter is |. As a corollary, we unconditionally obtain a non- 
trivial lower bound for the critical value of Hegarty 's random-orientation 
model. The second part of the paper is devoted to higher dimensions and 
we show that the Grimmett model percolates in any slab of height at least 
3 in Z 3 . 



1. Introduction 

Random-oriented percolation was first introduced by G. Grimmett HGrlll 
and is defined as follows. Consider the square lattice Z 2 and let each ver- 
tical edge be directed upwards with probability p e [0, 1] and downwards 
otherwise. Analogously, each horizontal edge is directed rightwards with 
probability p and leftwards otherwise. Let 6 G (p) be the probability that 
there is a directed path from the origin to infinity. By coupling with the 
classical bond percolation, it is not hard to show that 6g{\) — using, 
for example, Lemma 2.1 in [LinlJ. There is also the obvious symmetry 
9g(p) = #g(1 — p)> so it is natural to ask if 6 G (p) > for p ^ 1/2. This 
conjecture was raised by Grimmett HGrlB for the first time. The most signif- 
icant advance so far was also made by Grimmett in HGr2B where he showed 
that percolation does occur if one adds a positive density of randomly di- 
rected arcs, so that the total probability of a directed arc being present is 
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greater than one. Also, W. Xianyuan HXiall proved the uniqueness of the 
infinite cluster in the supercritical phase. He also conjectured that 9g(p) is 
strictly monotone on [1/2, 1], an obviously much stronger version of Grim- 
mett's conjecture. Both conjectures seem to be far from resolution. 

In this note we consider a slightly different model, proposed by P. Hegarty 
on MathOverflow [ |Hegl[ |, which we will refer as the //-model hereafter. It 
is defined as follows: for each edge e of the integer lattice Z 2 assign a 
direction away from the origin with probability p, and towards the origin 
otherwise. We say that a directed edge from x to y is oriented inwards if 
||x|| > \\y\\, outwards otherwise, with the usual Euclidean norm. We denote 
by Oh{p) the corresponding probability that there exists an infinite directed 
path from the origin. Note that this model coincides with the one proposed 
by Grimmett if we consider percolation only in the North-East quadrant. 
So, we will use 6 NE (p) for the percolation probability of the latter without 
abuse of notation. We conjecture that the NE quadrant is big enough for 
random-oriented percolation to occur. 

Conjecture 1. For random-oriented percolation in the North-East quad- 
rant, 9ne{p) > for all p > |. 

It might be possible to prove a weaker result that 6 G (p) > implies 
@ne(p) > 0, but we were unable to do so. The analogous result is well 
known for ordinary bond percolation in two dimensions and may be proven, 
for example, using RSW theory. 

It is not hard to show, using standard circuit-counting arguments, that 
0h(j>) = for p < -\, where // is the connective constant of the square 
lattice. From the other side, 9h(p) > for p > p c due to coupling with 
oriented percolation with the critical probability p c . It is proved that p c < 
0.6735, [BBS], and believed that p c w 0.6447. The main question that 
arises is whether we observe critical phenomena for the //-model and if so, 
what is the critical probability? Unfortunately, such a property appears to 
be very hard to establish. 

Conjecture 2. The probability function 6 nip) is strictly monotone in [0, 1]. 

In this note we prove 

Theorem 1. Suppose 6 NE {1 — p) > 0. Then 6 nip) = 0. 

Together with Conjecture [T| this would imply that for the //-model the 
critical probability does exist and it is equal to |. Also, we get the following 
result unconditionally. 

Corollary 1. 6 H {p) = for < p < 1 - p c . 
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Inserting the upper bound for oriented percolation, we get that 9h(p) = 
for p < 0.3265, which is considerably better than m 0.15. It is worth 
noting that the crucial property of the //-model is its' 90-degree rotational 
symmetry which is absent in the Grimmett model. 

At the end of the paper we consider the Grimmett model in higher di- 
mensions and prove that it is always supercritical, even if confined to a thin 
slab. 

Theorem 2. The 3-dimensional Grimmett model confined to the slab Z 2 x 
{ — 1, 0, 1} percolates for any p 6 [0, 1]. 

The proof of Theorem [2] exploits criticality of the two dimensional Grim- 
mett model that had already been shown in HGr2L Though the result sup- 
ports Grimmett's original conjecture, it seems that the crucial phenomena 
occur in the case of random-oriented percolation confined to a quadrant, 
which probably exhibits a phase transition in all dimensions. At least we 
can say that for any fixed d > 2 the //-model as well as the NE-quadrant 
model do not percolate for sufficiently small p > due to the standard 
path-counting argument, but of course they do percolate for p > p c . 

2. Percolation in the NE quadrant 

This section is devoted to the proof of Theorem [T] The dual lattice is 
a copy of Z 2 translated by the vector (1/2, 1/2), but orientation rules can 
be defined in two different ways: turning orientations in the original lattice 
clockwise or counterclockwise. We denote such dual lattices Z 2 / 1 and 7?/ 
and define them as follows. If the edge e fails to have an orientation in 
direction a, the dual edge has orientation a + n/2 in Z 2 / 1 and a — ti/2 in 
Z^ d . The corresponding dual NE quadrants we denote by and Qjjj. 

As an example, to prevent percolation in the NE quadrant, there must 
be a directed path from (x, —1/2) to (—1/2,?/) in and, equivalently, a 
directed path from (—1/2,?/) to (x, —1/2) in for some x,y > 0. That 
partly explains the superscripts u (up) and d (down). Also, we denote by 
A n the 2n x 2n square box with the center at the origin and let n (x) be 
the event that there exists a path from (sc, 1/2) to (1/2, y) in Q% for some 
y > which lies entirely within A m \ A n . For existence of a path which 
avoids A n we will write B+(x) and B + (x) for the unconstrained event. In 
other words, B+{x) = U m>n -B+ n (x) and B + (x) = U n> oB+(x). As usual, 
we will denote by dA n the vertex boundary of the box A n , i.e.: the set of 
vertices that have neighbors both inside and outside A n . 

Hereafter we will assume that all paths are in Z^ u if it is not explicitly 
stated otherwise. 

We start with a few auxiliary lemmas that explicitly exploit duality. 
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Lemma 1. If6 NE {l — p) > then 9ne(p) = 

Proof. This lemma can be proven in exactly the same way Harris showed 
there is no bond percolation in the quadrant at ~ in his seminal paper HHarrH . 
The only observation we need is that Z 2 / 1 has percolation parameter 1 — p 
if we fix its origin at some point (x, —1/2). Then, according to Lemma 5.2 
of HHarrlh since 6 NE {1 — p) > 0, with probability one there is an oriented 
path in from (x, —1/2) to (—1/2,?/) for some x,y > because any 
dual path started at the a>axis crosses the y-axis a.s. But this means that the 
NE-cluster in the original lattice is finite a.s. □ 

Lemma 2. Let n > and 6 = 9ne(1 — p) > 0. Recall that B+(x) denotes 
the event that there is a path in Q^from (x, 1/2) to (1/2, y) outside the box 
A n for some y > 0. Then 

liminfP{5+(x)} > 0. 

Proof. For each dual configuration tu, let u' be the modification of u> such 
that all edges inside A n are directed outwards from the point (1/2, 1/2). 
Let N(u) be the number of points (x, 1/2) such that tu E B + (x) but u' 
B + (x). Finally, let A be the set of all configurations tu, such that N(lo) = 
oo. We claim that F(A) = 0. 

Let us assume F(A) > for the sake of contradiction. For cu E A and x > 
0, conditions cu E B + (x) and cu' B + (x) imply existence of a path <9A n — >• 
(x — 1/2, y) in the original NE-quadrant of Z 2 for some y > 0. Indeed, 
since cu' ^ B + (x) there must be a NE-path in the original lattice that blocks 
(x, 1/2) from the y-axis in the dual (all other configurations would have 
probability zero). On the other hand, it can emanate only at the boundary 
of A n , because otherwise outwards orientation in A„ would have no effect 
on B + (x). But due to the fact that N(co) = oo we conclude that there 
must be arbitrarily long NE-paths in the original lattice, hence there exists 
an infinitely long one, implying 9ne(p) > F(A) > and contradicting 
Lemma [Q 

Now, defining N m (u) as above but counting only points (x, 1/2) with 
x > m we have 

lim F{{u\N m {u) > 0}) = 

m— >oo 

and thus 

liminf P{S+(x)} = liminf F{B + (x)} > 6. 

x— >oo x— >oo 

□ 

Corollary 2. Consider the H -model with edge probability p. Suppose 9 = 
^tvb(1— p) > 0. Then, for any n,d > OthereexistO < N < M < M such 
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that M — iV > d and 
for each x G [N, M ]. 

Proof. According to Lemma[2]we can pick iV such that F{B + (x) \ A n is blocked } > 
29/3 whenever x > N. As B+{x) = U m>n .B+ n (x) by definition, it re- 
mains to take Mo = N + d + 1 and M large enough to fulfill the desired 
conditions. □ 

Now we iterate Corollary E] to extend the directed path in the following 
way. Consider the event B° (A) that there exists a directed path 

(A, 1/2) (1/2, B) (-C, -1/2) (1/2, — D) -> (£, 1/2) (1) 

for some B,C, D, E > 0, where each part if the path, apart from axis 
crossings, lies inside a single dual quadrant. See Figure [TJ 

Lemma 3. For each N > there exists M > N such that 
F{B°(A) in A M \A N }> ^ ( ' l>V 



for some A e [N,M\. 

Proof. We may choose Mi > M > A^i > > such that F{B^ uN (x) > 
|} for all x G (JVi,M ). Thanks to Corollary S then, we pick M^,A^ 2 
(enlarging M and, subsequently, Mi if necessary) having M > Mq > 
A" 2 > Nx > such that P{S^ o 7Vi (z) > §} whenever x G {N 2 ,M' Q ). This 
guarantees that the probability of a directed path (x, 1/2) — >■ (1/2, B) — > 
(-C, 1/2) is greater than # 2 (1 — p)/4. Indeed, consider three events for any 
xG(iV 2 ,M^): 

(1) There exists a directed path (x, 1/2) — )■ (1/2, B) in Aa/ \ A^ 

(2) The axis crossing edge has direction (1/2, B) — > (—1/2, B) 

(3) There exists a directed path (-1/2,5) ->■ (-C, 1/2) inA Ml \ Ajy 
By construction, M > B > Ni and these three events are independent 
since they depend on disjoint edge sets. Thus, the probability that they oc- 
cur simultaneously is greater than 2 (1 —p)/A. Repeating this consideration 
twice, we get the claim of the corollary. Note, that here the rotational sym- 
metry of the //-model comes into play: turning by | each time we are able 
to construct the almost closed path with probability bounded away from 
zero. □ 

The following lemma asserts that the probability of a closed dual path 
is also bounded away from zero and is crucial. For A > 0, denote by 
B°(A, A) the event that there exists a closed directed circuit in Z^ u of the 
form ([[]) starting and finishing at (A, 1/2). We now claim 
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Lemma 4. For each N > there exists M > N such that 

F{B U (A, A) in A M \A N }> 1 - ^ ^ " 

for some A e [N,M\. 

Proof. Let us pick M, A > given by Lemma [3} such that 



F{B°(A) in A M \A N }> 



and fix the point A. Note that if there is at least one A — > B — > C — > 
D —7- E path, the inner-most and outer-most paths are unique and well- 
defined. Among all such paths (A, 1/2) -> (1/2,5) ->■ (-(7,-1/2) -> 
(1/2, —D) — > (E, 1/2) we choose the inner-most path P in and the outer- 
most path Pout, denoting their endpoints by E in and E out respectively, see 
Figure [T] Since the inner-most path always lies inside the outer-most one 
(though they may touch each other or even coincide), conditioning on the 
existence of least one A^B—^C—fD-^-E path, at least one of three 
events must occur: E in < A, E out > A or B° (A, A), whence 



F{E in <A} + F{E ont >A} + F(B°(A, A)) > 

If 

F(B U (A,A)) > 



i 



4 



3 V 2 

we are done, so we suppose 

PWo < A} > I 4 

without loss of generality. Conditioning on the states of the enclosed edges 
(the dashed region on the left part of Figure [T]), we ensure that the enclosing 
path is indeed inner-most. Now, recall the other dual lattice directed 
opposite to Z^ u . By symmetry, the probability for a clockwise oriented path 
(A, -1/2) (-1/2,-5') (-C',-1/2) -+ (—1/2, D') (P', 1/2) 
such that E' < A is equal to F{E- m < A] and any such path must cross 
(A, 1/2) -> (1/2, -> (-C in ,-l/2) -> (1/2, -An) -> (^,1/2) 
whenever E in < A. The states of the edges enclosed by the counter- 
clockwise path are independent of the existence of the clockwise path up 
to the first crossing point with the latter, and such a crossing produces a 
closed circuit in Z^ n (and, equivalently, in 1?/). Thus, 

F(B°(A, A)) > P{two paths exist, E in < A, E[ n < A} > ^ ~ P ^ . 
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Figure 1 . At least one of three events must occur: E- m < A 
(on the left), E ont > A (on the right) or B° (A, A) (closed 
circuit). In the first and second cases we condition on the 
states of edges in the dashed region to apply a symmetry 
argument. 



The case 

P[Emt >A} >i (fczr 4 

is completely analogous, but we condition on the states of the outer edges 
(see the right-hand side of Figure [T]). 

To finish the proof of Theorem [1} it simply remains to set up an infinite 
collection of frames provided by Lemma [4] and apply the Borell-Cantelli 
lemma. 

□ 



3. Higher dimensions 

In light of Theorem [T] it is reasonable to guess that both the H- and Grim- 
mett models are equivalent to percolation in the NE-quadrant. Indeed, sup- 
pose, say, 9 G (p) > and p > 0. According to Lemma [T] either 9ne(p) = 
or 9ne(1 ~ p) = (or maybe both). Analogously, due to self-duality of the 
NW and the SE quadrants of the Grimmett model for all p > it is easy to 
see that the percolation cluster restricted to any of these quadrants is almost 
surely finite. So, it is easy to believe that the infinite part of the percolation 
cluster (provided it is infinite) stays in the NE quadrant for p > 1/2, but we 
were unable to prove this. 

The (i-dimensional H- and Grimmett models differ significantly for d > 
3: whereas the if -model probably remains equivalent to percolation in the 
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quadrant, it is not difficult to show that 9q(p) > for all p when d > 3. In 
this section we will prove an even stronger result, namely that the Grimmett 
model percolates in any 3-dimensional slab of height at least three, as has 
already been announced in Theorem[2j On the other hand, the standard path 
counting argument implies that Q d H (p) = in all dimensions for sufficiently 
small p > (which, of course, depends on d). 

The idea of the proof of Theorem [2] is that one can consider spatial orien- 
tated paths of the form 

0, 2/, 0) -> (x, y, 1) ->■ (x + 1, y, 1) -> (x + 1, y, 0) 

as an additional arc (x, y, 0) — > (x+ l,y, 0) in the lattice Z 2 x {0} and then 
apply the following theorem by Grimmett HGr2B . 

Theorem 3. (Grimmett, HGr210 Consider the following independent pro- 
cess on Z 2 with parameters a and b: rightward and leftward ( respectively, 
upward and downward) arcs are placed independently between each pair 
of horizontal (respectively, vertical) neighbors. The probability of each up- 
ward or rightward arc being placed is a and the probability of each down- 
ward or leftward arc being placed is b. 

Ifa+b > 1 then the independent process with parameters a, b contains an 
infinite oriented self-avoiding path from with strictly positive probability. 

In the same paper HGr2H it is shown that the Grimmett model is equivalent 
to the independent process with parameters p and 1 — p, which implies that 
if additional arcs introduced above would have been placed independently, 
the process is supercritical. Thus, the main technical difficulty to overcome 
is the dependence between such paths for neighboring vertices in Z 2 x {0}. 

To do so, we will use Lemma 1.1 from ULSSH to "bound" the dependent 
measure by a product measure from below. 

Lemma 5. Suppose that (X s ) s( zs is a family of {0, l}-valued random vari- 
ables, indexed by a countable set S, with joint law v. Suppose S is to- 
tally ordered in such a way that, given any finite subset of S, si < s 2 < 
. . . < Sj < Sj + i, and any choice ofei, e 2 , . . . , ej G {0, 1}, then, whenever 
F(X Sl =e 1 ,...,X Sj =e j )>0, 

P(X Sj+1 = 1\X S1 = ei , . . . , X S] = ej ) > p. (2) 

Then v stochastically dominates tt^, which is a product measure with 
parameter p. 

Proof, (of Theorem [2]) We may assume that the original Bernoulli perco- 
lation has parameter p E [1/2, p c ] due to symmetry and obvious coupling 
with 2-dimensional directed percolation with the critical value p c . Chess- 
color the lattice vertices of Z 2 x {0} such that a vertex (x, y, 0) is black if 
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and only if x + y is even. For a given black vertex (x, y, 0) we consider 
oriented paths 

(x, y, 0) ->• (x, y, 1) ->• (x ± 1, y, 1) ->• (x ± 1, y, 0) (3) 

and 

(x, i/, 0) -> (x, y, 1) (x, y ± 1, 1) (x, y ± 1, 0). (4) 

For white vertices the construction is completely similar but all oriented 
paths go through the plane Z 2 x { — 1}. To make the distributions for white 
and black vertices identical and homogeneous we will consider a slightly 
different percolation model where each edge that is not in Z 2 x {0} has the 
same probability 1 — p < p of being oriented in any direction. This can be 
achieved by coupling in the following way: sample a countable set of in- 
dependent random variables uniformly distributed on [0, 1] for each edge in 
Z 2 x{ — 1,0,1} except the plane Z 2 x {0}. To make orientations distributed 
according to the original law we assign rightwards (resp. upwards) orienta- 
tion to the ith edge if Yi > 1 — p and leftwards (resp. downwards) otherwise. 
If we assign rightwards (resp. upwards) only if Yi E [2 — 2p, 1 — p) we will 
end up with the desired model with all orientations having the same proba- 
bility that is dominated by the original one. From now on we can claim that 
each auxiliary path is present with probability (1 — p) 3 . 

Let us fix the point (x, y, 0) for a moment and write Ai for the event that 
the i-th oriented path is present where the paths of the form ([3]) and ([4]) 
emanating from (x, y, 0) have been ordered in some way. We observe that 




for any (finite) index set /. On the other hand, each event Ai may be seen as 
an additional arc in Z 2 x {0} oriented outwards from (x, y, 0) and we end 
up with the probability measure v on Z 2 x {0} that corresponds to the set 
of arcs on Z 2 x {0} enriched in this way. 

Let us order the set of all possible additional arcs in S = Z 2 x {0} in 
some way, say, alphabetically. Note that, for each unoriented edge, there 
are two arcs directed in opposite ways and we consider them separately. 
Given the ordered countable set of arcs we assign a random variable X s . 
having X Si = 1 if the arc Si has been added during the enrichment process 
described above and zero otherwise. Obviously, the random variables X Si 
are not independent, but we are almost in the setting of Lemma [5] and it 
remains to show that inequality ([2]) holds for some p. 

Let us fix some X Si . First, we note that X s . is dependent only on arcs 
that either emanate from or end at the same vertex as Sj. Thus we have only 
six other variables X S1 , . . . , X s& upon which it depends. Moreover, due to 
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positive association ([5]) we have 

¥(X St = l\X n = ex, . . . , X 8 . = ej) > P(X Si = l\X n = 0, . . . , X S(i = 0) 

and it remains to bound the last probability from below. Without loss of 
generality we may assume that S{ is the arc (x, y, 0) — > (x + 1, y, 0). For 
both probabilities px,p 2 for arcs (x, y, 0) — > (x, y, 1) and (x + 1, y, 1) — >■ 
(x + l, 0) being present, conditionally on = 0, = 1, 6, we have a 
lower bound (1— p)p 3 because horizontal arcs lying on the planes Z 2 x {±1} 
are independent. Hence we can take p = (l—p) 3 p 6 , which is bounded away 
from zero since p G [1/2, p c ]. Thus, Lemma [5] applies and v dominates the 
product measure n p . 

But the process with the product measure n p corresponds to adding a 
positive density of additional arcs to Z 2 x {0} independently at random 
and we immediately arrive at the setting of Theorem |3j Thus, the original 
percolation process on Z 2 x { — 1, 0, 1} is supercritical. 

□ 

4. Conclusion 

Despite seemingly simple formulations, both Grimmett's and the i/-model 
appear to be difficult to analyze, mainly because of the absence of Harris- 
FKG-type correlation inequalities. The main difficulty to overcome is that 
any reasonable event, such as connectivity, defined in the random-orientation 
model is not increasing, in contrast with usual percolation models. Note, 
that Reimer's inequality still holds, but usually it is less fruitful and difficult 
to apply. It looks probable that further progress on the models in question 
requires substantially new ideas or at least considerable refinements of re- 
sults in classical percolation. 

In the proof of Theorem [T] we have shown how purely geometrical consid- 
erations based on rotational symmetry together with self-duality can be used 
even without joining paths by any kind of correlation inequalities. However, 
the conjecture that, for example, percolation in Grimmett's model implies 
percolation in the NE quadrant contains the Harris theorem (that there is 
no bond percolation at | in Z 2 ), and thus it is probably very non-trivial to 
prove. 

The crux of Theorem[2]consists of a few observations. First, as shown by 
Grimmett, the random-orientation process is equivalent to the independent 
process in which oriented arcs are placed independently. On the other hand, 
the process such that leftwards (resp. downwards) and rightwards (resp. 
upwards) orientations are independent and present with probabilities p and 
q respectively, dominates the same process with parameters p' < p and 
q' < q. Thus, it is monotone with respect to p and q and most of the classical 
results apply, for example, Menshikov's exponential decay theorem. 
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Second, in the same paper HGr2L Grimmett shows how self-duality (again 
without any correlation inequalities) together with exponential decay im- 
plies criticality (or, maybe, supercriticality) of the random-orientation model 
Finally, in the proof of Theorem [2] we use the general domination result to 
show that additional paths introduced by two copies of Z 2 , namely Z 2 x 
{ — 1} and Z 2 x {1}, have positive density and thus the resulting process is 
supercritical. Again, what we actually prove is that Grimmett's model in 
Z 2 x { — 1, 0, 1} dominates the independent process in Z 2 with some param- 
eters p' and q' such that p' + q' > 1. 
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